Efficiency of Information Spreading in a population of diffusing agents 
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We introduce a model for information spreading among a population of N agents diffusing on a 
square L x L lattice, starting from an informed agent (Source). Information passing from informed 
to unaware agents occurs whenever the relative distance is < 1. Numerical simulations show that the 
time required for the information to reach all agents scales as jV -a L' 3 , where a and /3 are noninteger. 
A decay factor z takes into account the degeneration of information as it passes from one agent to 
another; the final average degree of information of the population, I a v(z), is thus history-dependent. 
We find that the behavior of l av (z) is non-monotonic with respect to N and L and displays a set 
of minima. Part of the results are recovered with analytical approximations. 
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I. INTRODUCTION 

The information spreading in a population constitutes 
an attracting problem due to the emerging complex be- 
havior and to the great number of applications 
The propagation of information can be seen as a sequence 
of interpersonal processes between the interacting agents 
making up the system. In general, the population can 
be represented by a graph where agents are nodes and 
links between them exist whenever they interact with 
each other. 

Authors, who previously investigated the diffusion of 
information according to such a model, introduced differ- 
ent kinds of interpersonal interaction, but almost all of 
them assumed a static society 0, |H @ (a notable excep- 
tion being that of Eubank et al. In fact, networks 
are usually built according to a priori rules, which means 
that agents are fixed at their positions and can only in- 
teract with their (predetermined) set of neighbors (the 
flow of information between two agents is permanently 
open for linked pairs of agents and permanently closed 
for non- linked pairs). 

On the other hand, real systems are far from being 
static: nowadays individuals are really dynamic and con- 
tinuously come in contact, and lose contact, with other 
people. Hence, the interactions are rather instantaneous 
and time-dependent, and so should be considered the 
links of the pertaining graph. The network should be 
thought of as continuously evolving, adapting to the new 
interpersonal circumstances. 

Indeed, in sociology, where information spreading 
throughout a population is a long-standing problem 
it is widely accepted that processes of information trans- 
mission are far from deterministic. Rather, they should 
incorporate some stochastic elements arising, for exam- 
le, from "chance encounters with informed individuals" 
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Sociologists also underline that, irrespective of the kind 
of object to be transmitted, a realistic model should take 
into account whether the object passed from one agent 
to another is modified during the process [T^j. Espe- 



cially information, which spreads by replication rather 
than transference, is continuously revised while flowing 
throughout the network. Degradation during transmis- 
sion processes could reveal important qual itative and 
quantitative effects, as some recent works [Til fl2| j started 
to point out. 

This paper introduces a model that takes into account 
both the issues discussed above, namely, a mobile soci- 
ety and information changing during transmission. The 
model is based on a set of random walkers meant as "dif- 
fusing individuals" : a population of N interacting agents 
embedded on a finite space is represented by N random 
walkers diffusing on a square L x L lattice. We assume 
that two or more of them can interact if they are suffi- 
ciently close to each other: as a result, a given agent has 
no fixed position nor neighbors, but the set of agents it 
can interact with is updated at each instant. 

The information carried by an agent is a real (i.e., 
not boolean) variable, whose value lies between and 1. 
This (together with the diffusive dynamics) is the main 
point that differentiates our model from the susceptible- 
infected (SI) contact model of virus spreading in epidemi- 
ological literature 0], where only two status, susceptible 
and infected, are available to an agent. The issue of in- 
formation changing is dealt with by introducing a decay 
constant z < 1, which measures the corruption experi- 
enced by the piece of information when passing from an 
agent to another. We assume z to be universal: the more 
passages the information has undergone before reaching 
an individual, the more altered it is with respect to its 
original form. 

We study the time it takes for the piece of information 
to reach every agent (Population- Awareness Time). We 
show that it depends on N and L as a power-law, whose 
exponents are constant with respect to system parame- 
ters. We also investigate the final average (per agent) 
degree of information I av (z). We show that T av (z) is not 
a monotonic function of the density p = , but displays 
minima for definite values of N, L. This interesting re- 
sult implies that there does not exist a trivial direction 
where to tune the system parameters N and L in order 



2 



t = t = 1 t = 2 t = 3 




6 7 8 



FIG. 1: Evolution of 8 agents on a 5 x 5 lattice for t from to 3. For each t the lattice is shown above and the Information 
Tree is shown below. Informed agents are black circles; unaware agents are white circles. A grey circle of radius 1 is drawn 
around every informed agent to represent its action (an agent is in contact with another if it falls within this circle), t = 0: 
the only informed agent is the Information Source which carries information 1, so n(0, 0) = 1 and n(0, 1) = for / > 0. t = 1: 
agent 1 passes information to agent 2; now n(l,0) = 1, n(l, 1) = 1; t = 2: agent 1 passes information to agent 3 and agent 2 
passes information to agent 4; n(2, 0) = 1, n(2, 1) = 2, n(3, 1) = 1; t = 3: agent 2 passes information to agent 5; agent 4 passes 
information to agents 6, 7, 8. Notice that agent 6 is in contact with both 3 and 4; it chooses randomly to get information from 
4 (the same for agent 8). Now all agents have been informed: for this simulation the Population- Awareness Time is r = 3. The 
final information is T(r) — YlJ=o T ) zl = 1 + 2 z + 2 z 2 + 3 z 3 . 



to make information spreading more efficient. 

In the following, we first describe our model (Sec. HTj l. 
then we expose results obtained by means of numerical 
simulations (Sec. Illljl . Next, Sec. IIVI contains analyti- 
cal results which corroborate and highlight the former. 
Finally, Sec. is devoted to our conclusions and per- 
spectives. 



II. THE MODEL 

N random walkers (henceforth, agents) move on a 
square L x L lattice with periodic boundary condition. 
At time t = the agents are randomly distributed on 
the lattice. At each instant t > each agent jumps ran- 
domly to one of the four nearest-neighbor sites. There 
are no excluded- volume effects: there can be more agents 
on the same site; p = N/L 2 is the density of agents on 
the lattice. 

Each agent j carries a number Ij , Q < Ij < 1 , represent- 
ing information; an agent is called "informed" if Ij > 
and "unaware" if Ij = 0. At t = one agent, say agent 1, 
carries information 1 and is called the Information Source 
(or simply the Source); the other N — 1 agents are un- 
aware. The aim of the dynamics is to diffuse information 
from the Source to all agents. 

Interaction between two agents j and k takes place when 
i) one of them is informed and the other unaware, and ii) 
the chemical distance between the two agents is < 1 (i.e., 



they are either on the same site or on nearest-neighbor 
sites: we then say that they are "in contact"). By "in- 
teraction" we mean an information passing from the in- 
formed agent, say j, to the unaware one k with a fixed 
decay constant z (0 < z < 1): if j carries information Ij, 
then k becomes informed with information = z ■ 1%. 
Once an agent has become informed, it will never change 
nor lose its information (that is, informed agents never 
interact). If an unaware agent comes in contact with 
more informed agents at the same time, each carrying 
its own information Ij , it will acquire the information of 
one of them chosen at random (multiplied by z). The 
simulation stops at the time r when all the agents have 
become informed: we call this the Population- Awareness 
Time (PAT). 

We define n(t) the total number of informed agents at 
time t (n(0) = 1; n(r) = N). As a result of our model, 
the information carried by an agent is always a power of 
the decay constant z l , where I is the number of passages 
from the Information Source to the agent. We say that 
an informed agent belongs to level I when it has received 
information after I passages from the Information Source. 
We call n(l,t) the number of agents belonging to the l- 
th level at time t, or the population of level I at time t: 
n(t) = Ym=o n (l' t)- I n Fig- Qwe show as an example the 
evolution of = 8 agents on a 5 x 5 lattice. 
We can envisage information passing by drawing an In- 
formation Tree with N nodes and N — 1 links (fig. [IJ: 
the agents are the nodes of the tree, and a link is drawn 
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FIG. 2: (Color online) Evolution of n(t) for a population 
of N = 32 agents on six different lattices of size L — 2 m , 
m — 4, . . . , 9. Full circles denote the Population- Awareness 
Times, empty circles the Outbreak Times. 

between two agents when one passes information to the 
other. An agent belongs to level I if its distance from the 
Source along the tree is I. The Information Tree evolves 
with time: the tree at instant t is a subtree of that at 
instant t + 1. 

At each instant t we define the total information 

t 

I(z,f) = £n(I,t)z'; (1) 

1=0 

notice that it is the generating function of n(t); conse- 
quently, 

n(t) = 1(1, t). 

We are interested in particular in the final information 

J(z) =1(z,t). 
and in its average value per agent, I av (z) = I(z)/N. 

III. NUMERICAL RESULTS 

This section is divided into three parts. The first con- 
siders only n(t), the total informed population at time 
t, and the results presented are independent of the pop- 
ulation distribution on levels. The second section takes 
into account the distribution on levels n(l,t). The third 
section deals with the final information 2(z). All the re- 
sults are averaged over 500 different realizations of the 
system. 

A. Level-independent results 

Fig. [3 shows the typical time evolution of n(t), the 
number of aware people at time t, for fixed N and sev- 



eral different values of L. Due to the fact that, once 
informed, an agent can not modify his status, n(t) is a 
monotonic increasing function. The curve is sigmoidal: 
n(t) initially increases with an increasing growth rate 
dn(t)/dt. The growth rate is maximum at the Out- 
break Time tout, when usually n(t ou t) ~ N/2 (in Sec. 
IIVI we will justify this fact in a low-density approxima- 
tion). The growth rate then begins to decrease; the evo- 
lution slows down and the curve begins to saturate. The 
information reaches all the population at the Population- 
Awareness Time r, that is the quantity that we analyze 
here (roughly r ~ 2t ou t, and this fact as well will be 
justified in Sec. IIV() . 

The Population- Awareness Time r depends on the to- 
tal number of agents N and on the size of the lattice L, 
as shown in Fig. [31 As long as the density is not large 
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FIG. 3: (Color online) Dependence of the Population- 
Awareness Time r on the number of agents TV and the lattice 
size L. Top: Log- log scale plot of r versus iV; different lattice- 
size values are shown with different symbols and colors. For 
sufficiently small densities (p < 1), straight lines represent the 
best fit according to Eq. 01 Bottom: Log-log scale plot of r 
versus L; different values of the number of agents are shown 
with different symbols and colors. Provided that the density 
p is not large (p < 1), data points lay on the curves given by 
Eq. 21 which represent the best fit. Error on data points is 
<2%. 
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(p < 1)) data points are well fitted by power laws hold- 
ing over a wide range (though logarithmic corrections can 
not be ruled out): 

T~N~ a , (2) 
r~iA (3) 

The exponents a and /3 are constant by varying L or 
N, respectively, so that we can write: 

t <~ N~ a L@. (4) 

The fitting of data with an asymptotic least-squares 
method yields the following exponents: 

a = 0.68 ±0.01 (3 = 2.22 ±0.03. (5) 

B. Level-dependent results 

We now focus on the time evolution of n(l, t), the pop- 
ulation of level I. Each population evolves in time with a 
sigmoidal law (Fig. , with its own Outbreak Time and 
tending to a final value n(l, r). 

The final distribution of agents on levels n(l, r) as 
a function of I (Fig. 03 top) has an asymmetrical-bell 
shape, with a peak at position l pea k and a width <r, both 
depending on N and L (notice that only a fraction of the 
N available levels has a non- negligible population). If L 
is large enough (larger than L, see below), the population 
distribution on levels is well fitted by the 3-parameter 
function 

n(l,r) QogNf 
N T{B-l + C)' {> 

where T(x) is the Euler gamma function, and the param- 
eters A, B, C depend smoothly on N and L. The fitting 




levels 

FIG. 5: (Color online) Top: typical population distribution 
on levels at t = r for a low-density system (N=1024, L=4096). 
Squares are experimental results; the line is the result of data 
fitting according to Eq. ©. Bottom: population distribution 
on levels at t — t for systems with TV = 1024 and L between 
2 4 and 2 12 (the lines are guides to the eye): the behavior 
of the distribution is non-monotonic with respect to L. By 
increasing L from small values, the curves first shift to the 
right and flatten (L = 16, 24, 32). The rightmost, extremal 
curve corresponds to L — 64. Then, by increasing L the 
curves shift back to the left and sharpen (L = 128, 512, 4096). 




FIG. 4: (Color online) Time evolution of level populations 
n{l, t) for TV = 32, L = 512. 



function is a generalization of Eq. I|19|) , the distribution 
function of the low-density regime. 

In Fig. [SI bottom, we show how the distribution nil, r) 
changes with L for a fixed value N = 1024 and we in- 
troduce one of the most important results of this paper. 
For L small (hence for high density, p ^> 1) the distribu- 
tion is very sharp and peaked on small values of I. As 
L grows, the distribution shifts to higher values of I and 
becomes more and more spread (l pea k and a grow). The 
extremal, maximum-spread distribution is obtained for a 
value L — L (for N = 1024, L ~ 64): l pea k and a are 
maximum; the highest possible number of levels is occu- 
pied. As L is increased beyond L, the curve begins to 
shift back to smaller Is and to narrow; this process con- 
tinues up to the low-density regime (p <C 1). In general, 
L depends on N. 

The same phenomenon occurs if we keep L fixed and let 
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N vary. By increasing N from small, low-density values, 
the distribution shifts to the right and spreads, up to an 
extremal form occurring for N = N (depending on L). 
Then, it shifts back and narrows. 

This behavior has strong consequences on the efficiency 
of information spreading on the lattice, as we will see in 
the next section. 



C. Degree of Information 

In this section we deal with the final degree of informa- 
tion at the Population- Awareness Time, X{z) = T(r, z) 
(in particular, with its average value X av (z) — X(z)/N), 
and its dependence on N, L, and z. We remind (Eq. JJJ) 
that X{z) is the generating function of the final popula- 
tions n(l, t), hence its value depends on the final distri- 
bution of the population on levels analyzed in the previ- 
ous paragraphs. 

Once z is fixed, X av (z) depends nonmonotonically on N 
and L; let us follow it for N fixed and varying L in Fig.0 
For L small, due to the narrow distribution discussed 
in the previous section, the value of the information is 
high. When L = L, the population distribution on levels 
reaches its extremal form and the information displays a 
minimum. As L increases, the information starts to rise 
again. So, the main result is that, given a population 
number N, there is an optimal lattice size L for which 
the final information is minimum; this value is typically 
intermediate between the high-density and low-density 




FIG. 6: (Color online) Semilog scale plot of final degree of 
information per agent T av (z) = I(z)/N vs lattice size L. Sev- 
eral values of N are shown with different symbols and colors 
(lines are guides to the eye), while the decay constant is fixed 
at z = 0.9. Notice the occurrence of minima at L, N, and 




that L is monotonically increasing with respect to N. 
on data points is < 1.5%. 
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FIG. 7: (Color online) Semilog scale plot of final degree of 
information per agent I av (z) as a function of the lattice size 
L, when N = 2 9 (lines are guides to the eye). Four different 
values of decay constant z are considered, as shown by the 
legend. 



regimes. The same happens having fixed L and letting TV 
vary: there is a minimum for N — N, where N depends 
on L. 

This result implies that choosing an optimization strat- 
egy for the spreading of information on the lattice is not 
trivial. Suppose e.g. that we are given N agents on a 
lattice and we want to maximize the final average infor- 
mation X av {z) by varying the lattice size L (starting from 
some Lq). This optimization process is meant to be local: 
we are not allowed to modify the size by several orders of 
magnitude, but just around the starting size Lq. Then, 
the choice whether to shrink or expand the lattice de- 
pends on Lq. If Lq < L, increasing L takes the system 
closer to the information minimum (X av (z) decreases); 
decreasing L increases X av (z) and is the right strategy. 
If on the other hand Lq > L, increasing L is the right 
strategy. 

Fig. \7\ shows that the depth of the information mini- 
mum depends in turn on the decay constant z: as z is 
varied from to 1, there are some curves (corresponding 
to in-between values) which display a more emphasized 
minimum. 

Finally, in Fig.|S|we show how the final average degree 
of information X av (z) depends on z, for different values 
of N, once the size L is fixed. There are, as expected, 
two fixed points: when z = 1 [z = 0), X av (z) is equal to 1 
(0), irrespective of the parameters (N,L) of the system. 
The function X av (z) cannot be determined but in two 
particular regimes (low- and high-density). 

When p = N/L 2 is sufficiently low (p < 2~ 8 ), the 
function is well fitted by 

X av (z) = N z -\ (7) 

within the error (< 3%). When p > 1, X av (z) is fitted by 
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FIG. 8: (Color online) Final (t — r) degree of information 
per agent T av (z) versus the decay constant z. The size of 
the lattice is fixed as L — 2 4 , while several values of JV are 
considered and represented in different colors and symbols. 
The curve depicted is the best fit when N = 2 10 (p > 1), 
according to Eq. |S] Notice the existence of the fixed points 
2 = 0, T av (z) = jj and z = l,X av (z) = 1. 




L-1 levels 

FIG. 9: (Color online) Evolution of the system in the high- 
density approximation for a lattice with L = 12 and a Source 
starting in the site with coordinates (7, 7). Left: the wave 
front of information on the lattice at times t = 1, 2, 3, 4 
has a square shape. Agents in the interior of the square are 
informed, agents on the exterior are unaware. Times corre- 
spond to levels: agents between the front at time 3 and that 
at time 4 belong to level 4, and so on. Right: final distribu- 
tion of agents on levels; broken lines highlight the triangular 
shape of the distribution. 



(1 - z B - L ) 2 

l av {z)=A-zp±— (8) 

(1 - zf 

with A, B depending on N, L. 

The two laws come from particular population distri- 
butions, as will be explained in the next section. 

IV. ANALYTICAL RESULTS 

Consider a system with N and L fixed. Let V(t) be the 
probability that at time t an unaware agent is in contact 
with at least 1 informed agent. Let Vi(k,s; t) be the 
probability that at time t an unaware agent is in contact 
with k + s informed agents, of which k belonging to level 
I and s belonging to some other level. Then the evolution 
of the system is governed by two master equations, one 
for the total population: 

n(t + l) =n(t) + (N-n(t))V{t), (9) 

and one for the level populations: 

n(l, t + 1) = n(l, t) + (N- n{t)) V Vi-i{k, s; t)-^—. 

t— 1 k + s 

k,s 

(10) 

V(t) and Vi{k, s; t) are very complex functions of their 
arguments and cannot be calculated in the general case. 
For example, V(t) depends not only on the number of in- 
formed agents n(t) but also on their spatial distribution, 
hence on the instant and the site where each of them 
has been informed (in other words, on the history of the 



system). We will calculate the evolution of the system 
in two particular cases, for high and low densities, and 
finally compare the results with intermediate systems. 

High-density regime. In this case (p — ► oo) there are 
many agents on every site. If the agents on a site get in- 
formed at a time t, we can suppose that at t + 1 at least 
one of them will jump on each of the four nearest- neigh- 
bor sites: hence, all the unaware agents on the nearest- 
and next-to-nearest-neighbor sites will get information 
at time t + 1. In this way (Fig. ^ information spread- 
ing among agents amounts to propagation of information 
through the lattice. A "wave front" of information trav- 
els with constant velocity: on the interior sites are in- 
formed agents, on the exterior sites unaware agents. If 
we suppose the Source to be at the center of the lat- 
tice at t = 0, at each instant the wave front is the lo- 
cus of points whose chemical distance from the center 
is 2t + 1. Consequently, n(t) — p(8t 2 — 4t + 1), up to 
the half- filling time t ou t ~ when the front reaches 
the boundary of the lattice; for t > L/4, the equation is 
n(t) = p(~8t 2 + U{2L + 1) + (L + l) 2 ). The Population- 
Awareness Time is r ~ L/2. 

Almost all the agents on the wave front at time t have 
received information at time t — 1: so, each new time 
step adds a new level, whose population never changes at 
successive times. The population n(l, t) is proportional 
to the length of the wave front at the time t — I: n(l, t) ~ 
4p(4Z + 1) up to t = L/A and n(l, t) - 4p(-4l + 2L - 1) 
up to t = L/2. As can be seen from Fig. the shape 
of the level distribution at t — r is triangular (compare 
this to the distribution for L = 16 in Fig. [SJ. The Final 
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Information is proportional to p, according to the formula 

N L/4 T 

I(z) =^2n(l, t)z 1 ~^16pZz' + Yl 4p(2i-4Z). 

1=0 1=0 l=L/4+l 

(1 _ z L/4 )2 



16zp- 



(11) 



A modified version of this equation, Eq.©, has 
been used to fit the information curves for high-density 
regimes. 

Low-density regime. In the case of low density (p <C 1) 
the time an informed agent walks before meeting an un- 
aware agent becomes very large. We can then assume 
that the agents between each event have the time to re- 
distribute randomly on the lattice, that is, we adopt a 
mean- field approximation. Let p — 5/L 2 be the proba- 
bility that two given agents, randomly positioned on the 
lattice, are in contact (5 is the number of points contained 
in a circle of radius 1). Hence, (1 —p) n w is the probabil- 
ity for an agent at time t of not being in contact with any 
of the n(t) informed agents, and V(t) = 1 — (1 — p) n ^> 
is the probability of being in contact with at least one 
informed agent. Master equation ij^J) becomes: 

n(t + 1) = n(t) + (N- n{t)) (l - (1 - p)"W) , 



and to first order in p: 



n(t + 1) = n(t) +p {N - n(t)) n(t). 



(12) 



Thus, n(t + 1) = f(n(t)); f is a logistic-like map, with 
a repelling fixed point in (/'(0) = 1 + Np), and an 
attracting fixed point in N (f'(N) = 1 — Np). Since 
Np — 5p <C 1, the increment of n(t) at each time step is 
very small (of order p), and we can take the evolution to 
be continuous. The equation becomes: 



n(t+ 1) -n(t) 



dn(t) 
dt 



p(N -n(t))n(t) (13) 



and the solution, with the initial condition n(0) = 1, is 
the sigmoidal function 



i(t) = N 



a Npt 



e Npt + N 



(14) 



The outbreak time, i.e. the flex of the curve, is in t out = 

log(iV-l) 



Np 



that is also the half-filling time, n(t c 



N/2. 



The total population N is reached only for t — oo, but 
we can take the PAT to be the time when N — 1 agents 
have been informed: 



21og(JV-l) 21ogN 



Np 



Np 



fogiV 2 

N 



(15) 



where the last result holds for N large: hence, in the 
low-density approximation the exponent for L is f3 — 2, 
while the law for N contains logarithmic corrections and 



the exponent a cannot be defined. The first result in Eq. 
(|15fl shows that in this approximation r = 2 t out ■ 
The quantity Vi{k,s; t) in Eq. (UJl is: 

*(M;t)= (»(J:*))(n(*)-»(M) 
xp k+s (l - 

The sum over k and s in Eq. l|l(Jf) . using the Chu- 
Vandermonde identity for binomial coefficients [T^| . 
yields a master equation for the level populations in the 
mean-field approximation: 

n(l, t+1) - n(l, t ) + (N-n(t))(l-(l- P r^) n{l ~ *' * } , 

n[t) 

and to first order in p: 

n(l, t + 1) =n(l, t)+pn(l-l, t) (N - n(t)). 

Its continuous version is: 
dn(l, t) 



dt 



pn{l-l, t)(N-n(t)) 



(16) 



that has to be solved for each I. For 1 = 1, with the initial 
condition n(l, 0) = 0, we get the solution 

n(l,t) = Npt- \og(e Npt + N- l) + logN 
= log Hi)). 

We then plug this solution into Eq. (|16fl to get n(2, t), 
and so on. It can be shown by induction that for every 
I, with the initial condition n(/, t) = 0, 

n{l, t) =-{Npt- log (e Npt + N - l) + logA)' 

= ^(n(l,t)) l = ±[log(n(t))] 1 . (17) 



This set of curves (not shown here) is similar to that of 
Fig. 01 with crossovers and different Outbreak Times. 
The normalized level population at each t is: 

n(U)_ 1 1 l = e^^n m l 



n(t) n(t) 11 



l\ 



( 18 ) 

hence, it is a Poisson distribution with mean log (n(t)). 
The population distribution on levels at t — r is 



n(l, t) 



l\ 



(19) 



independent of p (hence of L). A modified version of this 
distribution, Eq. ©, has been used to fit the numerical 
curves. 

The total information is 

N N 

l(t, z) =J2 n(h t)z l = £ - [log (n(tj) ■ zf ~ 



~ e 



i=0 1=0 
\og(n(t))-z 



n(ty. 



(20) 
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L= 16 L= 32 




L=64 L=512 

FIG. 10: Snapshots of four systems with N = 1024 and 
L — 16, 32, 64, 512, all at an instant near to the half-filling 
time. Only informed agents are shown; they are represented 
as circles of radius 1. Notice that the high-density picture of 
a connected set of informed agents holds up to L = L — 64. 
For L > 64, the picture breaks down and the system is better 
described by a low-density approximation (L = 512). 

In particular, T(t, z) = N z , in agreement with Eq. J7J). 

In conclusion, we have examined the system in two 
different regimes, both optimal for information spread- 
ing. The worst case for information spreading, at L, 



seems to correspond to crossover between these two 
regimes, as shown in Fig. I1UI 



V. CONCLUSIONS AND PERSPECTIVES 

We have presented a model of information spreading 
amongst diffusing agents. The model takes into account a 
population made up of agents who are socially, as well as 
geographically, dynamic. Moreover, it allows for possible 
alteration of information occurring during the transmis- 
sion process, by introducing a decay constant z. 

Investigations are lead both by means of numerical 
simulations and of analytical methods valid in the high- 
and low-density regimes. 

The main results are two. First: the time it takes the 
piece of information to reach the whole population of N 
agents, distributed on a lattice sized L, depends on N and 
L according to a power law. This behavior holds over a 
wide range, where exponents are found to be constant 
and noninteger. Second: the final (i = r) average degree 
of information X av (z) for a fixed population N (lattice 
size L) shows a surprisingly non-monotonic dependence 
on the lattice size L (on the population N), with the 
occurrence of a minimum. This means that, from an 
applied perspective, an optimization strategy for T av (z) 
is possible with respect to N and L. 

Extensions of our model to networks embedded in 
topologically different spaces are under study. 
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We introduce a model for information spreading among a population of N agents diffusing on a 
square L x L lattice, starting from an informed agent (Source). Information passing from informed 
to unaware agents occurs whenever the relative distance is < 1. Numerical simulations show that the 
time required for the information to reach all agents scales as jV -a L' 3 , where a and /3 are noninteger. 
A decay factor z takes into account the degeneration of information as it passes from one agent to 
another; the final average degree of information of the population, I a v(z), is thus history-dependent. 
We find that the behavior of l av (z) is non-monotonic with respect to N and L and displays a set 
of minima. Part of the results are recovered with analytical approximations. 
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I. INTRODUCTION 

The information spreading in a population constitutes 
an attracting problem due to the emerging complex be- 
havior and to the great number of applications 
The propagation of information can be seen as a sequence 
of interpersonal processes between the interacting agents 
making up the system. In general, the population can 
be represented by a graph where agents are nodes and 
links between them exist whenever they interact with 
each other. 

Authors, who previously investigated the diffusion of 
information according to such a model, introduced differ- 
ent kinds of interpersonal interaction, but almost all of 
them assumed a static society 0, |H @ (a notable excep- 
tion being that of Eubank et al. In fact, networks 
are usually built according to a priori rules, which means 
that agents are fixed at their positions and can only in- 
teract with their (predetermined) set of neighbors (the 
flow of information between two agents is permanently 
open for linked pairs of agents and permanently closed 
for non- linked pairs). 

On the other hand, real systems are far from being 
static: nowadays individuals are really dynamic and con- 
tinuously come in contact, and lose contact, with other 
people. Hence, the interactions are rather instantaneous 
and time-dependent, and so should be considered the 
links of the pertaining graph. The network should be 
thought of as continuously evolving, adapting to the new 
interpersonal circumstances. 

Indeed, in sociology, where information spreading 
throughout a population is a long-standing problem 
it is widely accepted that processes of information trans- 
mission are far from deterministic. Rather, they should 
incorporate some stochastic elements arising, for exam- 
le, from "chance encounters with informed individuals" 
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Sociologists also underline that, irrespective of the kind 
of object to be transmitted, a realistic model should take 
into account whether the object passed from one agent 
to another is modified during the process [T^j. Espe- 



cially information, which spreads by replication rather 
than transference, is continuously revised while flowing 
throughout the network. Degradation during transmis- 
sion processes could reveal important qual itative and 
quantitative effects, as some recent works [Til fl2| j started 
to point out. 

This paper introduces a model that takes into account 
both the issues discussed above, namely, a mobile soci- 
ety and information changing during transmission. The 
model is based on a set of random walkers meant as "dif- 
fusing individuals" : a population of N interacting agents 
embedded on a finite space is represented by N random 
walkers diffusing on a square L x L lattice. We assume 
that two or more of them can interact if they are suffi- 
ciently close to each other: as a result, a given agent has 
no fixed position nor neighbors, but the set of agents it 
can interact with is updated at each instant. 

The information carried by an agent is a real (i.e., 
not boolean) variable, whose value lies between and 1. 
This (together with the diffusive dynamics) is the main 
point that differentiates our model from the susceptible- 
infected (SI) contact model of virus spreading in epidemi- 
ological literature 0], where only two status, susceptible 
and infected, are available to an agent. The issue of in- 
formation changing is dealt with by introducing a decay 
constant z < 1, which measures the corruption experi- 
enced by the piece of information when passing from an 
agent to another. We assume z to be universal: the more 
passages the information has undergone before reaching 
an individual, the more altered it is with respect to its 
original form. 

We study the time it takes for the piece of information 
to reach every agent (Population- Awareness Time). We 
show that it depends on N and L as a power-law, whose 
exponents are constant with respect to system parame- 
ters. We also investigate the final average (per agent) 
degree of information I av (z). We show that T av (z) is not 
a monotonic function of the density p = , but displays 
minima for definite values of N, L. This interesting re- 
sult implies that there does not exist a trivial direction 
where to tune the system parameters N and L in order 
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FIG. 1: Evolution of 8 agents on a 5 x 5 lattice for t from to 3. For each t the lattice is shown above and the Information 
Tree is shown below. Informed agents are black circles; unaware agents are white circles. A grey circle of radius 1 is drawn 
around every informed agent to represent its action (an agent is in contact with another if it falls within this circle), t = 0: 
the only informed agent is the Information Source which carries information 1, so n(0, 0) = 1 and n(0, 1) = for / > 0. t = 1: 
agent 1 passes information to agent 2; now n(l,0) = 1, n(l, 1) = 1; t = 2: agent 1 passes information to agent 3 and agent 2 
passes information to agent 4; n(2, 0) = 1, n(2, 1) = 2, n(3, 1) = 1; t = 3: agent 2 passes information to agent 5; agent 4 passes 
information to agents 6, 7, 8. Notice that agent 6 is in contact with both 3 and 4; it chooses randomly to get information from 
4 (the same for agent 8). Now all agents have been informed: for this simulation the Population- Awareness Time is r = 3. The 
final information is T(r) — YlJ=o T ) zl = 1 + 2 z + 2 z 2 + 3 z 3 . 



to make information spreading more efficient. 

In the following, we first describe our model (Sec. HTj l. 
then we expose results obtained by means of numerical 
simulations (Sec. Illljl . Next, Sec. IIVI contains analyti- 
cal results which corroborate and highlight the former. 
Finally, Sec. is devoted to our conclusions and per- 
spectives. 



II. THE MODEL 

N random walkers (henceforth, agents) move on a 
square L x L lattice with periodic boundary condition. 
At time t = the agents are randomly distributed on 
the lattice. At each instant t > each agent jumps ran- 
domly to one of the four nearest-neighbor sites. There 
are no excluded- volume effects: there can be more agents 
on the same site; p = N/L 2 is the density of agents on 
the lattice. 

Each agent j carries a number Ij , Q < Ij < 1 , represent- 
ing information; an agent is called "informed" if Ij > 
and "unaware" if Ij = 0. At t = one agent, say agent 1, 
carries information 1 and is called the Information Source 
(or simply the Source); the other N — 1 agents are un- 
aware. The aim of the dynamics is to diffuse information 
from the Source to all agents. 

Interaction between two agents j and k takes place when 
i) one of them is informed and the other unaware, and ii) 
the chemical distance between the two agents is < 1 (i.e., 



they are either on the same site or on nearest-neighbor 
sites: we then say that they are "in contact"). By "in- 
teraction" we mean an information passing from the in- 
formed agent, say j, to the unaware one k with a fixed 
decay constant z (0 < z < 1): if j carries information Ij, 
then k becomes informed with information = z ■ 1%. 
Once an agent has become informed, it will never change 
nor lose its information (that is, informed agents never 
interact). If an unaware agent comes in contact with 
more informed agents at the same time, each carrying 
its own information Ij , it will acquire the information of 
one of them chosen at random (multiplied by z). The 
simulation stops at the time r when all the agents have 
become informed: we call this the Population- Awareness 
Time (PAT). 

We define n(t) the total number of informed agents at 
time t (n(0) = 1; n(r) = N). As a result of our model, 
the information carried by an agent is always a power of 
the decay constant z l , where I is the number of passages 
from the Information Source to the agent. We say that 
an informed agent belongs to level I when it has received 
information after I passages from the Information Source. 
We call n(l,t) the number of agents belonging to the l- 
th level at time t, or the population of level I at time t: 
n(t) = Ym=o n (l' t)- I n Fig- Qwe show as an example the 
evolution of = 8 agents on a 5 x 5 lattice. 
We can envisage information passing by drawing an In- 
formation Tree with N nodes and N — 1 links (fig. [IJ: 
the agents are the nodes of the tree, and a link is drawn 
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FIG. 2: (Color online) Evolution of n(t) for a population 
of N = 32 agents on six different lattices of size L — 2 m , 
m — 4, . . . , 9. Full circles denote the Population- Awareness 
Times, empty circles the Outbreak Times. 

between two agents when one passes information to the 
other. An agent belongs to level I if its distance from the 
Source along the tree is I. The Information Tree evolves 
with time: the tree at instant t is a subtree of that at 
instant t + 1. 

At each instant t we define the total information 

t 

I(z,f) = £n(I,t)z'; (1) 

1=0 

notice that it is the generating function of n(t); conse- 
quently, 

n(t) = 1(1, t). 

We are interested in particular in the final information 

J(z) =1(z,t). 
and in its average value per agent, I av (z) = I(z)/N. 

III. NUMERICAL RESULTS 

This section is divided into three parts. The first con- 
siders only n(t), the total informed population at time 
t, and the results presented are independent of the pop- 
ulation distribution on levels. The second section takes 
into account the distribution on levels n(l,t). The third 
section deals with the final information 2(z). All the re- 
sults are averaged over 500 different realizations of the 
system. 

A. Level-independent results 

Fig. [3 shows the typical time evolution of n(t), the 
number of aware people at time t, for fixed N and sev- 



eral different values of L. Due to the fact that, once 
informed, an agent can not modify his status, n(t) is a 
monotonic increasing function. The curve is sigmoidal: 
n(t) initially increases with an increasing growth rate 
dn(t)/dt. The growth rate is maximum at the Out- 
break Time tout, when usually n(t ou t) ~ N/2 (in Sec. 
IIVI we will justify this fact in a low-density approxima- 
tion). The growth rate then begins to decrease; the evo- 
lution slows down and the curve begins to saturate. The 
information reaches all the population at the Population- 
Awareness Time r, that is the quantity that we analyze 
here (roughly r ~ 2t ou t, and this fact as well will be 
justified in Sec. IIV() . 

The Population- Awareness Time r depends on the to- 
tal number of agents N and on the size of the lattice L, 
as shown in Fig. [31 As long as the density is not large 
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FIG. 3: (Color online) Dependence of the Population- 
Awareness Time r on the number of agents TV and the lattice 
size L. Top: Log- log scale plot of r versus iV; different lattice- 
size values are shown with different symbols and colors. For 
sufficiently small densities (p < 1), straight lines represent the 
best fit according to Eq. 01 Bottom: Log-log scale plot of r 
versus L; different values of the number of agents are shown 
with different symbols and colors. Provided that the density 
p is not large (p < 1), data points lay on the curves given by 
Eq. 21 which represent the best fit. Error on data points is 
<2%. 
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(p < 1)) data points are well fitted by power laws hold- 
ing over a wide range (though logarithmic corrections can 
not be ruled out): 

T~N~ a , (2) 
r~iA (3) 

The exponents a and /3 are constant by varying L or 
N, respectively, so that we can write: 

t <~ N~ a L@. (4) 

The fitting of data with an asymptotic least-squares 
method yields the following exponents: 

a = 0.68 ±0.01 (3 = 2.22 ±0.03. (5) 

B. Level-dependent results 

We now focus on the time evolution of n(l, t), the pop- 
ulation of level I. Each population evolves in time with a 
sigmoidal law (Fig. , with its own Outbreak Time and 
tending to a final value n(l, r). 

The final distribution of agents on levels n(l, r) as 
a function of I (Fig. 03 top) has an asymmetrical-bell 
shape, with a peak at position l pea k and a width <r, both 
depending on N and L (notice that only a fraction of the 
N available levels has a non- negligible population). If L 
is large enough (larger than L, see below), the population 
distribution on levels is well fitted by the 3-parameter 
function 

n(l,r) QogNf 
N T{B-l + C)' {> 

where T(x) is the Euler gamma function, and the param- 
eters A, B, C depend smoothly on N and L. The fitting 




levels 

FIG. 5: (Color online) Top: typical population distribution 
on levels at t = r for a low-density system (N=1024, L=4096). 
Squares are experimental results; the line is the result of data 
fitting according to Eq. ©. Bottom: population distribution 
on levels at t — t for systems with TV = 1024 and L between 
2 4 and 2 12 (the lines are guides to the eye): the behavior 
of the distribution is non-monotonic with respect to L. By 
increasing L from small values, the curves first shift to the 
right and flatten (L = 16, 24, 32). The rightmost, extremal 
curve corresponds to L — 64. Then, by increasing L the 
curves shift back to the left and sharpen (L = 128, 512, 4096). 




FIG. 4: (Color online) Time evolution of level populations 
n{l, t) for TV = 32, L = 512. 



function is a generalization of Eq. I|19|) , the distribution 
function of the low-density regime. 

In Fig. [SI bottom, we show how the distribution nil, r) 
changes with L for a fixed value N = 1024 and we in- 
troduce one of the most important results of this paper. 
For L small (hence for high density, p ^> 1) the distribu- 
tion is very sharp and peaked on small values of I. As 
L grows, the distribution shifts to higher values of I and 
becomes more and more spread (l pea k and a grow). The 
extremal, maximum-spread distribution is obtained for a 
value L — L (for N = 1024, L ~ 64): l pea k and a are 
maximum; the highest possible number of levels is occu- 
pied. As L is increased beyond L, the curve begins to 
shift back to smaller Is and to narrow; this process con- 
tinues up to the low-density regime (p <C 1). In general, 
L depends on N. 

The same phenomenon occurs if we keep L fixed and let 
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N vary. By increasing N from small, low-density values, 
the distribution shifts to the right and spreads, up to an 
extremal form occurring for N = N (depending on L). 
Then, it shifts back and narrows. 

This behavior has strong consequences on the efficiency 
of information spreading on the lattice, as we will see in 
the next section. 



C. Degree of Information 

In this section we deal with the final degree of informa- 
tion at the Population- Awareness Time, X{z) = T(r, z) 
(in particular, with its average value X av (z) — X(z)/N), 
and its dependence on N, L, and z. We remind (Eq. JJJ) 
that X{z) is the generating function of the final popula- 
tions n(l, t), hence its value depends on the final distri- 
bution of the population on levels analyzed in the previ- 
ous paragraphs. 

Once z is fixed, X av (z) depends nonmonotonically on N 
and L; let us follow it for N fixed and varying L in Fig.0 
For L small, due to the narrow distribution discussed 
in the previous section, the value of the information is 
high. When L = L, the population distribution on levels 
reaches its extremal form and the information displays a 
minimum. As L increases, the information starts to rise 
again. So, the main result is that, given a population 
number N, there is an optimal lattice size L for which 
the final information is minimum; this value is typically 
intermediate between the high-density and low-density 




FIG. 6: (Color online) Semilog scale plot of final degree of 
information per agent T av (z) = I(z)/N vs lattice size L. Sev- 
eral values of N are shown with different symbols and colors 
(lines are guides to the eye), while the decay constant is fixed 
at z = 0.9. Notice the occurrence of minima at L, N, and 




that L is monotonically increasing with respect to N. 
on data points is < 1.5%. 



Error 



FIG. 7: (Color online) Semilog scale plot of final degree of 
information per agent I av (z) as a function of the lattice size 
L, when N = 2 9 (lines are guides to the eye). Four different 
values of decay constant z are considered, as shown by the 
legend. 



regimes. The same happens having fixed L and letting TV 
vary: there is a minimum for N — N, where N depends 
on L. 

This result implies that choosing an optimization strat- 
egy for the spreading of information on the lattice is not 
trivial. Suppose e.g. that we are given N agents on a 
lattice and we want to maximize the final average infor- 
mation X av {z) by varying the lattice size L (starting from 
some Lq). This optimization process is meant to be local: 
we are not allowed to modify the size by several orders of 
magnitude, but just around the starting size Lq. Then, 
the choice whether to shrink or expand the lattice de- 
pends on Lq. If Lq < L, increasing L takes the system 
closer to the information minimum (X av (z) decreases); 
decreasing L increases X av (z) and is the right strategy. 
If on the other hand Lq > L, increasing L is the right 
strategy. 

Fig. \7\ shows that the depth of the information mini- 
mum depends in turn on the decay constant z: as z is 
varied from to 1, there are some curves (corresponding 
to in-between values) which display a more emphasized 
minimum. 

Finally, in Fig.|S|we show how the final average degree 
of information X av (z) depends on z, for different values 
of N, once the size L is fixed. There are, as expected, 
two fixed points: when z = 1 [z = 0), X av (z) is equal to 1 
(0), irrespective of the parameters (N,L) of the system. 
The function X av (z) cannot be determined but in two 
particular regimes (low- and high-density). 

When p = N/L 2 is sufficiently low (p < 2~ 8 ), the 
function is well fitted by 

X av (z) = N z -\ (7) 

within the error (< 3%). When p > 1, X av (z) is fitted by 
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FIG. 8: (Color online) Final (t — r) degree of information 
per agent T av (z) versus the decay constant z. The size of 
the lattice is fixed as L — 2 4 , while several values of JV are 
considered and represented in different colors and symbols. 
The curve depicted is the best fit when N = 2 10 (p > 1), 
according to Eq. |S] Notice the existence of the fixed points 
2 = 0, T av (z) = jj and z = l,X av (z) = 1. 




L-1 levels 

FIG. 9: (Color online) Evolution of the system in the high- 
density approximation for a lattice with L = 12 and a Source 
starting in the site with coordinates (7, 7). Left: the wave 
front of information on the lattice at times t = 1, 2, 3, 4 
has a square shape. Agents in the interior of the square are 
informed, agents on the exterior are unaware. Times corre- 
spond to levels: agents between the front at time 3 and that 
at time 4 belong to level 4, and so on. Right: final distribu- 
tion of agents on levels; broken lines highlight the triangular 
shape of the distribution. 



(1 - z B - L ) 2 

l av {z)=A-zp±— (8) 

(1 - zf 

with A, B depending on N, L. 

The two laws come from particular population distri- 
butions, as will be explained in the next section. 

IV. ANALYTICAL RESULTS 

Consider a system with N and L fixed. Let V(t) be the 
probability that at time t an unaware agent is in contact 
with at least 1 informed agent. Let Vi(k,s; t) be the 
probability that at time t an unaware agent is in contact 
with k + s informed agents, of which k belonging to level 
I and s belonging to some other level. Then the evolution 
of the system is governed by two master equations, one 
for the total population: 

n(t + l) =n(t) + (N-n(t))V{t), (9) 

and one for the level populations: 

n(l, t + 1) = n(l, t) + (N- n{t)) V Vi-i{k, s; t)-^—. 

t— 1 k + s 

k,s 

(10) 

V(t) and Vi{k, s; t) are very complex functions of their 
arguments and cannot be calculated in the general case. 
For example, V(t) depends not only on the number of in- 
formed agents n(t) but also on their spatial distribution, 
hence on the instant and the site where each of them 
has been informed (in other words, on the history of the 



system). We will calculate the evolution of the system 
in two particular cases, for high and low densities, and 
finally compare the results with intermediate systems. 

High-density regime. In this case (p — ► oo) there are 
many agents on every site. If the agents on a site get in- 
formed at a time t, we can suppose that at t + 1 at least 
one of them will jump on each of the four nearest- neigh- 
bor sites: hence, all the unaware agents on the nearest- 
and next-to-nearest-neighbor sites will get information 
at time t + 1. In this way (Fig. ^ information spread- 
ing among agents amounts to propagation of information 
through the lattice. A "wave front" of information trav- 
els with constant velocity: on the interior sites are in- 
formed agents, on the exterior sites unaware agents. If 
we suppose the Source to be at the center of the lat- 
tice at t = 0, at each instant the wave front is the lo- 
cus of points whose chemical distance from the center 
is 2t + 1. Consequently, n(t) — p(8t 2 — 4t + 1), up to 
the half- filling time t ou t ~ when the front reaches 
the boundary of the lattice; for t > L/4, the equation is 
n(t) = p(~8t 2 + U{2L + 1) + (L + l) 2 ). The Population- 
Awareness Time is r ~ L/2. 

Almost all the agents on the wave front at time t have 
received information at time t — 1: so, each new time 
step adds a new level, whose population never changes at 
successive times. The population n(l, t) is proportional 
to the length of the wave front at the time t — I: n(l, t) ~ 
4p(4Z + 1) up to t = L/A and n(l, t) - 4p(-4l + 2L - 1) 
up to t = L/2. As can be seen from Fig. the shape 
of the level distribution at t — r is triangular (compare 
this to the distribution for L = 16 in Fig. [SJ. The Final 
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Information is proportional to p, according to the formula 

N L/4 T 

I(z) =^2n(l, t)z 1 ~^16pZz' + Yl 4p(2i-4Z). 

1=0 1=0 l=L/4+l 

(1 _ z L/4 )2 



16zp- 



(11) 



A modified version of this equation, Eq.©, has 
been used to fit the information curves for high-density 
regimes. 

Low-density regime. In the case of low density (p <C 1) 
the time an informed agent walks before meeting an un- 
aware agent becomes very large. We can then assume 
that the agents between each event have the time to re- 
distribute randomly on the lattice, that is, we adopt a 
mean- field approximation. Let p — 5/L 2 be the proba- 
bility that two given agents, randomly positioned on the 
lattice, are in contact (5 is the number of points contained 
in a circle of radius 1). Hence, (1 —p) n w is the probabil- 
ity for an agent at time t of not being in contact with any 
of the n(t) informed agents, and V(t) = 1 — (1 — p) n ^> 
is the probability of being in contact with at least one 
informed agent. Master equation ij^J) becomes: 

n(t + 1) = n(t) + (N- n{t)) (l - (1 - p)"W) , 



and to first order in p: 



n(t + 1) = n(t) +p {N - n(t)) n(t). 



(12) 



Thus, n(t + 1) = f(n(t)); f is a logistic-like map, with 
a repelling fixed point in (/'(0) = 1 + Np), and an 
attracting fixed point in N (f'(N) = 1 — Np). Since 
Np — 5p <C 1, the increment of n(t) at each time step is 
very small (of order p), and we can take the evolution to 
be continuous. The equation becomes: 



n(t+ 1) -n(t) 



dn(t) 
dt 



p(N -n(t))n(t) (13) 



and the solution, with the initial condition n(0) = 1, is 
the sigmoidal function 



i(t) = N 



a Npt 



e Npt + N 



(14) 



The outbreak time, i.e. the flex of the curve, is in t out = 

log(iV-l) 



Np 



that is also the half-filling time, n(t c 



N/2. 



The total population N is reached only for t — oo, but 
we can take the PAT to be the time when N — 1 agents 
have been informed: 



21og(JV-l) 21ogN 



Np 



Np 



fogiV 2 

N 



(15) 



where the last result holds for N large: hence, in the 
low-density approximation the exponent for L is f3 — 2, 
while the law for N contains logarithmic corrections and 



the exponent a cannot be defined. The first result in Eq. 
(|15fl shows that in this approximation r = 2 t out ■ 
The quantity Vi{k,s; t) in Eq. (UJl is: 

*(M;t)= (»(J:*))(n(*)-»(M) 
xp k+s (l - 

The sum over k and s in Eq. l|l(Jf) . using the Chu- 
Vandermonde identity for binomial coefficients [T^| . 
yields a master equation for the level populations in the 
mean-field approximation: 

n(l, t+1) - n(l, t ) + (N-n(t))(l-(l- P r^) n{l ~ *' * } , 

n[t) 

and to first order in p: 

n(l, t + 1) =n(l, t)+pn(l-l, t) (N - n(t)). 

Its continuous version is: 
dn(l, t) 



dt 



pn{l-l, t)(N-n(t)) 



(16) 



that has to be solved for each I. For 1 = 1, with the initial 
condition n(l, 0) = 0, we get the solution 

n(l,t) = Npt- \og(e Npt + N- l) + logN 
= log Hi)). 

We then plug this solution into Eq. (|16fl to get n(2, t), 
and so on. It can be shown by induction that for every 
I, with the initial condition n(/, t) = 0, 

n{l, t) =-{Npt- log (e Npt + N - l) + logA)' 

= ^(n(l,t)) l = ±[log(n(t))] 1 . (17) 



This set of curves (not shown here) is similar to that of 
Fig. 01 with crossovers and different Outbreak Times. 
The normalized level population at each t is: 

n(U)_ 1 1 l = e^^n m l 



n(t) n(t) 11 



l\ 



( 18 ) 

hence, it is a Poisson distribution with mean log (n(t)). 
The population distribution on levels at t — r is 



n(l, t) 



l\ 



(19) 



independent of p (hence of L). A modified version of this 
distribution, Eq. ©, has been used to fit the numerical 
curves. 

The total information is 

N N 

l(t, z) =J2 n(h t)z l = £ - [log (n(tj) ■ zf ~ 



~ e 



i=0 1=0 
\og(n(t))-z 



n(ty. 



(20) 
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L= 16 L= 32 




L=64 L=512 

FIG. 10: Snapshots of four systems with N = 1024 and 
L — 16, 32, 64, 512, all at an instant near to the half-filling 
time. Only informed agents are shown; they are represented 
as circles of radius 1. Notice that the high-density picture of 
a connected set of informed agents holds up to L = L — 64. 
For L > 64, the picture breaks down and the system is better 
described by a low-density approximation (L = 512). 

In particular, T(t, z) = N z , in agreement with Eq. J7J). 

In conclusion, we have examined the system in two 
different regimes, both optimal for information spread- 
ing. The worst case for information spreading, at L, 



seems to correspond to crossover between these two 
regimes, as shown in Fig. I1UI 



V. CONCLUSIONS AND PERSPECTIVES 

We have presented a model of information spreading 
amongst diffusing agents. The model takes into account a 
population made up of agents who are socially, as well as 
geographically, dynamic. Moreover, it allows for possible 
alteration of information occurring during the transmis- 
sion process, by introducing a decay constant z. 

Investigations are lead both by means of numerical 
simulations and of analytical methods valid in the high- 
and low-density regimes. 

The main results are two. First: the time it takes the 
piece of information to reach the whole population of N 
agents, distributed on a lattice sized L, depends on N and 
L according to a power law. This behavior holds over a 
wide range, where exponents are found to be constant 
and noninteger. Second: the final (i = r) average degree 
of information X av (z) for a fixed population N (lattice 
size L) shows a surprisingly non-monotonic dependence 
on the lattice size L (on the population N), with the 
occurrence of a minimum. This means that, from an 
applied perspective, an optimization strategy for T av (z) 
is possible with respect to N and L. 

Extensions of our model to networks embedded in 
topologically different spaces are under study. 
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